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(ii) the simulation procedure of the dechanneling process of a positron in periodically 
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of the undulator over broad ranges of the photons. The calculations are performed for 
various parameters of the channel bending. 
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1. Introduction 

In this paper we proceed with the investigation of the properties of the spectrum of 
emitted photons accompanying ultra-relativistic positron planar channeling through a 
crystal which is periodically bent, as it is illustrated in figure [TJ This phenomenon was 
described recently in ^ |2| and was called Acoustically Induced Radiation (AIR). It 
was noted that the periodic pattern of crystal bendings (which can be achieved either 
through propagation of a transverse acoustic wave or by using static periodically strained 
crystalline structures @, §]) gives rise to a new mechanism of electromagnetic emission 
of the undulator type, in addition to a well-known ordinary channelling radiation [|], 
The mechanism of the AIR is as follows. Provided certain conditions are fulfilled |2| the 
beam of particles, which enters the crystal at a small incident angle with respect to the 
curved crystallographic plane, penetrates through the crystal following the bendings of 
its channel. This results in transverse oscillations of the beam particles (additional to 
the oscillations inside the channel due to the action of the interplanar force). These 
oscillations become an effective source of spontaneous radiation of undulator type 
due to the constructive interference of the photons emitted from similar parts of the 
trajectory. It was demonstrated that the system "ultra-relativistic charged particle 
+ periodically bent crystal" serves as a new type of undulator, and, consequently, as a 
new source of undulator radiation of high intensity, monochromaticity and of a particular 
pattern of the angular-frequency distribution. 

The main subject of the present study is the detailed and accurate quantitative 
consideration of the influence of the dechanneling process [BJ on the parameters of the 
AIR which are the intensity and the specific pattern of the spectral-angular and spectral 
distributions. The dechanneling, i.e. the decrease in the volume density of the channeled 
particles n(z) with penetration distance z due to the multiple scattering of the projectiles 
with the target electrons and nuclei, is the parasitic effect which leads to the restriction 
on the crystal length L and, correspondingly, on the number of the undulator periods 
N, which, in turn, defines the intensity of the AIR radiation |fj]. The factor, which 
plays a crucial role in obtaining the accurate data for the AIR characteristics, is a 
so-called dechanneling length Ld, which is a mean penetration distance covered by a 
channeling particle. The quantity La depends on the projectile energy, the type of the 
crystal and crystallographic plane, and on the parameters of the channel bending (these 
are explained in more detail below in this section). In this connection the following 
problems are solved and discussed below in the paper: 

(i) We propose the scheme for accurate quantitative treatment of the AIR in presence 
of the dechanneling (section |3|). As a result, we evaluate simple analytic expression 
for spectral-angular distribution of the AIR which contains, as a parameter, the 
dechanneling length. 
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(ii) We introduce the algorithm of the simulation procedure of the dechanneling process 
of a positron in periodically bent crystals (section |]). 

(iii) Using (ii) we calculate the dechanneling lengths of 5 GeV positrons channeling in 
Si, Ge and W crystals the (110) plane of which are periodically bent (section |5]). 
The calculations are performed for various parameters of the channel bending. 

(iv) Using (i) and (iii) we calculate the spectral-angular and spectral distributions of 
the AIR over broad ranges of the photons emitted by 5 GeV positrons channeled 
in Si and W (section []). 
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Figure 1. Schematic representation of the periodically bent crystallographic plane. 
The open circles mark the atoms belonging to two neighbouring crystallographic planes 
(which are parallel to the (xz) plane) in the initially linear crystal (d is the interplanar 
spacing). The filled circles denote the positions of the atoms when the crystal is 
periodically bent. The z axis is directed along the straight channel center line, the y 
axis marks the direction perpendicular to the crystallographic planes in the straight 
crystal. The profile of the bent channel centerline (the dashed curve) is described by a 
periodic function y = S(z) the amplitude, a, and the period, A, of which satisfy a«A. 



In our previous papers (see also [|7|, ||) we have estimated the role of the 
dechanneling by assuming simple relation L < to be valid. The magnitude of 
Ld was estimated on the basis of the diffusion theory (see ||, [TC^, 11, 12]) with some 
corrections introduced due to the bending of the crystal |L2L [L3| . On the basis of these 



estimates we carried out the calculations of the AIR spectra for various energies of 
the projectile, different types of the crystals and crystallographic plane and for various 
parameters of the channel bending (these are explained in more detail below). It was also 
pointed out |2| that the system "ultra-relativistic charged particles + periodically bent 
crystal" leads, in addition to the spontaneous radiation, to the possibility of generating 
stimulated emission, similar to the one known for a free-electron laser jnj in which the 
periodicity of a trajectory of an ultra-relativistic projectile is achieved by applying a 



4 



spatially periodic magnetic field. In connection with the stimulated AIR it was noted 
in |2| that to achieve a noticeable degree of amplification one has to operate with a 
positron bunch of a high volume density. The dechanneling length played a crucial role 
in these estimations by influencing the intensity of the spontaneous AIR as well as the 
gain factor and the corresponding value of the positron bunch volume density for the 
stimulated AIR. Indeed, increasing of the dechanneling length by a factor of 2 results 
in decreasing of the channeling beam density needed to achieve large gain almost by by 
a factor of 10 0. 

By utilizing the numeric scheme outlined in section f| (to our knowledge it is the 
first time that the dechanneling process is simulated in a periodically bent crystal) we 
demonstrate further in section ||] that the simple estimation of Ld made in J| in many 
cases produces the results which are noticeably smaller then more accurate ones. This 
results, consequently, in essential (up to the order of magnitude for heavy crystals) 
difference in the intensities of the AIR radiation (section ||) computed by using the 
estimated, L e d , and the calculated, L d , values of the dechanneling lengths. 

Prior to presenting our main results let us formulate the important criteria which 
are used further in the paper. 

The planar channeling of an ultra-relativistic projectile in a straight crystal (for 
reviews see |15|, [16], [IT]) occurs only if the entrance angle of the projectile velocity with 



respect to the mid-plane does not exceed the critical angle of channeling || 

e < e L (i) 

where the quantity Ql (often called as the Lindhard angle) is defined as (2U max /m'y). 
Here U max is the interplanar potential-well depth, 7 is the relativistic factor of the 
projectile and m is its mass (below we consider only the channeling of positrons, if 
otherwise is not explicitly stated). 

The criterion for a stable channeling of an ultra-relativistic particle in a bent crystal 



was formulated in ]18[ and has clear physical meaning: the maximum centrifugal force 
due to the channel bending must be less than that of the interplanar field. Provided 
this condition is fulfilled, the beam of channeling particles at each instant moves inside 
the channel, especially parallel to the bent crystal midplane as it does in the case 
of linear channeling. The dynamics of the particles channeled in crystals bent with a 
constant (or slowly varied) curvature 1/R is described in details in [^, |lj, [H], ^ 

(see also the references therein). In the cited reviews the main accent was made on 
channeling of heavy projectiles due to the growing interest for using bent crystals for 
extraction, bending and splitting of high-energy beams (for the latest review see 
and the references therein. 

For a periodically bent crystal the criterion for a stable channeling can be written 
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in the form || 

0= 5 -^-<l, (2) 

-"-mm ^max 

where ?7 max is the maximum value of the interplanar field, -R m i n is the minimum curvature 
radius of the bent channel centerline. Introducing the shape function S(z), which 
describes the profile of the centerline (see figure |1]), the minimum curvature l/_R mm can 
be written in terms of the amplitude a and the period A of S(z): 1/-R mm ~ \S"(z)\^ x ~ 
ci~ 1 (\/2tt) 2 . Thus, for a given crystal and its crystallographic plane, the condition 2 
establishes the ranges of projectile energies and the parameters a and A inside which 
the channeling in a periodically bent channel can effectively occur. 

The channeling radiation |4], [| formed by ultra-relativistic positrons channeled 
in various straight crystals has been studied extensively both theoretically (e.g. 



17] , 16 , 23), and references therein) and experimentally [[25|, [27], |28|, |29fl . 

Several papers |J0[ [H], [2TJ] have treated the radiation emitted by light projectiles 
channeled in non-periodically one-arc bent crystals. In this case the channel bending 
gives rise to the synchrotron-type radiation in addition to the channeling one. In the 
cited papers it was demonstrated that for small values of the parameter C (i.e. for 
large curvatures of the bending) the photon energies corresponding to the maximum 
intensity of the synchrotron radiation do not exceed the energy of the first harmonic of 
the channeling radiation. 

When the shape of the bent channel becomes periodic the additional radiation 



becomes of the undulator type |32|, [HI, E3, R3[ rather than the synchrotron one. Then 



there appear two essentially different regimes which are defined by the magnitude of 
the ratio a/d |,|. In the case a/d <C 1, which can be realized either by applying 
low-amplitude transverse ultrasonic wave |3~|, B5|, |37|, or by using a supperlattice 
crystalline structure the characteristic frequencies of the two types of radiation can 



become compatible, and thus, one can consider, at least theoretically, the phenomenon of 
the resonant enhancement of the photon yield due to the coupling of these mechanisms. 
However, as it was demonstrated in || |7], H , * n ^ ne nm h a/d <C 1 the intensity of AIR is 
noticeably small compared not only with that of the channeling radiation but also with 
the intensities of the radiative background mechanisms, such as coherent and incoherent 
bremsstrahlung. Hence, it is highly questionable whether the AIR can be considered as 
a new phenomenon in the limit of low values of a. 

On the contrary, if the amplitude of the shape function and the interplanar spacing 
satisfy strong inequality 

a > d (3) 

then (i) the characteristic frequencies of the AIR and the channeling radiation are well- 
separated (this is the consequence of both conditions, (3) and C < 1), (ii) the intensity 
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of the AIR is essentially higher that that of the ordinary channeling radiation |2], [7], |8|]. 
In this case the undulator-type radiation due to the periodic structure of the crystal 
bending can be considered as a new source of the emission within the X- and 7-range. 

Below in the paper we assume that both strong inequalities (2) and (3) are fulfilled. 
It is important to mention that in our scheme the amplitude a is subject not only to 
the condition (3) but also is much smaller compared with the period of S(z): a < A. 
Hence, neither the interplanar spacing nor the distance between neighbouring lattice 
atoms are changed (in more detail it is discussed in section pTT| ). 

To conclude the introductory part we note that two realistic ways of "preparation" of 
a periodically bent channel may be discussed |l|, @]. ^ ^ s feasible, by means of modern 
technology, to grow the crystal with its channels been statically bent according to a 
particular pattern. In particular, the crystalline undulator can be constructed based 
on graded composition strained layers in a superlattice || (see also [^, |40|). Another 
possibility is to use high-amplitude transverse ultrasonic wave. In the latter case the 
parameters of the undulator can be easily tuned by varying the amplitude and the 
frequency of the acoustic wave. It was demonstrated that the acoustically-based 
undulator can be created by applying the supersonic waves within the frequency range 
10 . . . 100 MHz which is achievable in the experiments on propagating the positron beams 



through acoustically-excited crystals EI 



2. The AIR radiation in absence of the dechanneling 

Let us consider the spectral-angular distribution of the photons emitted during the 
projectile positron planar channeling in a periodically bent crystal along the z direction 
(see figure 0). The motion of the particle occurs in the (yz)-plane, so that its radius 
vector and the velocity are given by r(t) = ze z + ye y , v(t) = ze z + ye y where e z and 
e y are the unit vectors along the z- and y-directions, respectively. 

Provided the conditions (|l]-3) are fulfilled the photon spectrum is formed due to 
the ordinary channeling mechanism and to the AIR, which prevail over two other 
radiative background mechanisms, coherent and incoherent bremsstrahlung. It has 



been established, both theoretically and experimentally (e.g. [16|, |28|), that the 
total radiative losses of a positron channeling in a straight crystal are essentially smaller 
than its energy provided e < 10 — 20 GeV. This estimate is correct in the case of 
channeling through periodically bent crystal as well 0. 

Therefore, restricting from above the range of e by the value ~ lOGeV, one can 
assume that the emitted photon energies satisfy hu <C e. This inequality allows 
to consider the process of the photon emission within the framework of classical 
electrodynamics (e.g. 



The distribution of the energy radiated in the cone df2 n along the direction n and 
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within the frequency interval [u, u + da;] is given by the following expression 

4 = , 4((**^« (4) 

where a = e 2 /hc ~ 1/137, r is the time of flight through the crystal of thickness L, 
t m L/ c. The solid angle f2 n of the photon emission is characterized by two polar angles 
•& and ip. The functions 4>{ti,t 2 ) and f(t 1 ,t 2 ) are given by 

<f>(t u t 2 )=t 1 -t 2 --n-(r 1 -T 2 ), f{h,t 2 ) = - (p^p.-l^ (5) 

The short-hand notations used here are r 1 )2 = r(fi )2 ) and v 12 = v(t 12 ). 

Assuming that the relativistic factor of the projectile satisfies strong inequality 
7 ^> 1, one can carry out the expansion of the functions <f>(ti, t 2 ) and f(ti, t 2 ) in powers 
of 7 _1 (see e.g. ||16|| ). Then, retaining the senior non- vanishing terms, one transforms 
the r.h.s. of (H) to a form which is more suitable for analytical and numerical analysis: 

_^ = fia ^{|/ 1 -.co S ,/ o p + # sm V|/o| 2 } (6) 



where 



The modified phase function <p(t) reads as 



dte^, I^Tdt^e^ (7) 



, 1 ^ 2 p 2 \ A(t) 

^)=* (^ + Y + ^)+4 i -^ co ^^ U (8) 



A(t)= dW^-^ , (9) 




P 2 = ^ % (10) 

The last relation defines the so-called undulator parameter p, which is commonly used 
in the theory of the undulator radiation (e.g. |16|). This quantity is related to the mean- 
square velocity y 2 in the ^/-direction, which is transverse with respect to the centerline 
of the bent channel (see figure 0). 

If one neglects the processes of inelastic scattering by the crystal electrons and 
nuclei, then the motion of the positron in the channeling regime can be described within 
the framework of the continuum approximation || [Tj|] for the interaction potential U 



between the projectile and lattice atoms. Within this approximation and in the case 
of periodically bent channel the dependence y(t), which enters the r.h.s. of (0-9) and, 
thus, defines the spectral-angular distribution (|6|), can be obtained as follows 0, ||. 

At any given moment t the ^/-coordinate of the projectile can be presented as a sum 
of two terms 

y(t) = S(z)+p(t) (11) 
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Here S(z) is the displacement due to the channel bending. The quantity p(t) stands for 
the displacement of the projectile from the centerline of the channel. The equation of 
motion for p(t) reads 

p = -—-c 2 S" (12) 

where the dot sign stands for the time derivative, U' = dU / dp is the derivative of the 
interplanar potential which depends on p, and S" = d 2 S"/d2; 2 . The first term on the 
r.h.s. represents the acceleration of the particle due to the action of the interplanar 
force, while the second one is the acceleration due to the centrifugal force. Note that 
on the r.h.s. one may put z = ct after carrying out the differentiation S", because all 
the deviations from c of the particle's velocity along the z-axis were taken into account 
when writing eqs. (0-9). 

Solving the equation of motion (12) and, then, using the result, combined with 
(11), in (0-9), one obtains the total spectral- angular distribution of the radiation @. 
This scheme was used in our recent papers pi where the results of the numerical 
calculations of the spectra were presented for various shapes S(z) and with using 
different approximations for the interplanar potential U(p). 

A particle which channels through periodically bent crystal experiences two types 
of oscillations in the transverse direction, as it is seen from (11). Firstly, there are the 
channeling oscillations due to the action of the interplanar potential. The oscillations 
of the second type (henceforth we refer to them as to the "undulator oscillations") are 
related to the periodicity in the shape S(z) of the centerline of the channel inside which 
the particle moves. The characteristic frequencies of these oscillations, denoted as Q c h 
and fl u respectively, can be estimated as follows 



/ U' 2nc 

fidi ~ i/3 j = — 7- (13) 

V arwy A 

Correspondingly, the quantities Q c h, Q u define the characteristic frequencies, u c h and 
u u , of the photons emitted due to the channeling and the undulator oscillations: 

uj ch ~ 7 2 tt ch , u u ~ 7 2 tt u (14) 

Taking into account the relations (2) and (3) one estimates the relative magnitude of 
Q ch and Q u : 

ii ch 10 ch art mm^ m ax a 

Hence, the characteristic frequencies of the AIR and the ordinary channelling radiation 
are well separated J/J, and, what is also very important, in the region u ~ u u there is 
no coupling of the two mechanisms of the radiation H. 
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The results of numerical calculation of the spectral distribution (|6|) presented in || 
clearly suggest that the characteristics of the AIR spectra (the position of the peaks 
of the radiation, their width, the radiated intensity) are, practically, insensitive to the 
choice of the approximation used to describe the interplanar potential U(p). Moreover, 
if being interested in the spectral- angular distribution only in the region u ~ u u , then 
one may disregard the channeling oscillations and to assume that the projectile moves 
along the centerline of the bent channel. In this case, taking into account the periodicity 
of the undulator motion, the r.h.s. of (Q) can be presented in the form which clearly 
exhibits all the characteristic features of the AIR radiation. 

Let assume that there is integer number of the undulator periods on the crystal 
length, L = A N. Also, for the sake of transparency of the analytic expressions, let 
consider the shape function S(z) of the form 



S(z) = a sinkz, k = 2n/\ 



(16) 



Let assume that there is integer number of the undulator periods on the crystal 
length, L = XN. Then, the spectral-angular distribution of the AIR radiation can 
be represented in the form [EL E4| 



dE 



N 



du dQ n 



u 2 4tt 2 
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1? COS p Fq 



+ d 2 sin 2 v? |F | 



(17) 



Here the subscript N in the notation dE^ indicate the number of undulator periods, 
uq = 2nc/X is the so-called undulator frequency, and the quantities F m {m = 0, 1) stand 
for the integrals 



V 




x 



r/ijj + — sin(2-0) — (3 sin ip 



x 



uj p 
u>o 47 2 



/3 = i? cos ip 



uj p 
uj 7 



(19) 



The information on the specific pattern of the undulator radiation, namely on the 
positions of the narrow maxima in the spectral-angular distribution (|17D , is accumulated 
in the factor D^{rf) which equals to 



Dn(v) 



( sin Ntttj^ 
\ sin ni] 



(20) 



This function, which is well-known in the theory of undulator radiation (see e.g. |T6|, |24|) 
and in classical diffraction theory [^2] , defines the profile of the line of the emission in an 
ideal undulator. It has sharp maxima at the points 77 = k — 1, 2, 3... where D^(k) = N 2 . 
The integer values rj = k define (see the first equation in (19)) the characteristic 
frequencies Uk (harmonics) of the undulator radiation for which the spectral-angular 
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distribution flTTT) reaches its maxima: 



47 2 to k 

Uk = 2 + 2^ +p 2 (21) 

The natural width of each line of the emission can be estimated as Au/ = (l/N)(u>k/k) 
and is independent on the harmonic number k. 

The largest number of the radiated harmonics, k max , can be estimated as k max ~ p 3 /2 
(see e.g. 0, ) , so that in the range to < co , / Cmax the spectral-angular distribution of the 
AIR represents by itself the set of narrow well-separated peaks. The radiation intensity 
at u) = u k is proportional to the square of the total number of the undulator periods, 
N 2 , reflecting the coherence effect of radiation. 

It follows from fl2"T|) that the position of the peaks are dependent on the relativistic 
factor, on the emission angle and on the undulator parameter p. 

The only feature, related to the channeling oscillations, which modifies the formulae 
(0-21), which describe the AIR spectrum, concerns the definition of the undulator 
parameter p (see (10)). 

When the channeling oscillations are disregarded completely (this means omitting 
the term p(t) on the r.h.s. of (11)) then the undulator parameter equals 2 / -f 2 S' 2 = 
47r 2 7 2 a 2 /A 2 , and for given 7 is defined only by the amplitude and the period of the 
shape function S(z). 

To account for the channeling oscillations and to modify the definition of p one can 
use the following arguments. The quantity p 2 is proportional to the the mean-square 
velocity y 2 . The latter, according to (11), can be written as follows: 

p = c 2 ^ 2 + ~p 2 + 2cW~p 

The last term on the r.h.s. contains the product of the factors oscillating with 
incompatible frequencies (see (|1~5D, and, therefore, it disappears when one carries out the 
averaging over the time interval ~ 1/Q c h- Hence, the square of the undulator parameter 
can be written as follows 

P 2 =l 2 (^) 2 + 7 2 ^ (22) 

The results presented in [[| demonstrate that in the region uj ~ uj u the spectral-angular 
distribution, obtained by means of the exact procedure which includes (i) numerical 
integration of the equation of motion (12) with the interplanar potential considered 
within the the Moliere approximation, (ii) subsequent evaluation of the integrals in 
(0), and (iii) by further averaging of the result over all channeling trajectories, can be 
accurately reproduced by using (19) and (22) with the term p 2 calculated within the 
harmonic approximation for U(p). The use of (19) instead of simplifies the numerical 
procedures considerably, leading to the reduction, by orders of magnitude, of the CPU 
time. 
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3. The AIR radiation in presence of the dechanneling 

If the dechanneling is neglected, one may unrestrictedly increase the intensity of the AIR 
by considering larger N- values. In reality, random scattering of the channeling positron 
by the electrons and nuclei of the crystal leads to a gradual increase of the particle energy 
associated with the transverse oscillations in the channel. As a result, the transverse 
energy at some distance from the entrance point exceeds the depth of the interplanar 
potential well, and the particle leaves the channel. Consequently, the volume density 
n(z) of the channeled particles decreases with the penetration distance z. Although the 
exact explicit dependence of the channeled fraction of the particles n(z)/no on z (no is 
the volume density of the channeled particles at the entrance) hardly can be obtained 
by analytical means due to the complexity of the accurate treatment of the multiple 
scattering problem in medium, it has been argued 0, [l7j that far from the entrance 
point n(z)/n can be described by the exponential decay law 

= exp (-z/L d ) (23) 

n 

Basing on this relation, below in this section we present a model which allows to 
modify the spectral-angular distribution (|17D by taking into account the influence of 
the dechanneling effect. 

Let us assume that the period A of the shape function S(z) and the dechanneling 
length L d satisfy the strong inequality A <C L d . If so, then the number of particles (per 
unit volume) An = n(z) — n(z+\), which dechannel within the interval [z, z+X], is small 
compared with the number of the particles n(z) which reach the point z. The inequality 
allows also to introduce the (discrete) probability pj of the event, that after channeling 
through j periods of the undulator (j = 1, 2, ... N — 1) the particle dechannels within 
the subsequent (j + l)-th period. In accordance with (|23|) the probability pj is written 
as follows 

= f exp (-j/N d ) (exp (1/N d ) - 1) , for j = 1, 2 ... N - 1 
P ° \exp(-(N-l)/N d ) for j = N 

and is subject to the normalization condition J2jL\Pj = 1- in (|24l) we introduced 
the quantity N d = L d /X which is the number of the undulator periods within the 
dechanneling length. Since N d 3> 1, then the deviation of N d from the nearest integer is 
of no importance, and henceforce the quantity N d is treated as the integer much larger 
than one. 

With the help of (|24|) one can construct the following quantity, which characterizes 
the spectral-angular distribution of the AIR radiation formed in the crystal of the length 
L (and of the corresponding number of the undulator periods N) and which accounts 
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for the dechanneling effect: 

/ dE N \ = 
\ da; df2 n / 



N 



da; df2 n 

(^(77)) 



4vr 2 



|- 

7 



$ cosy?F | 2 + $ 2 sinV|F | 2 



(25) 
(26) 



with 



N 



(exp (l/iV d ) - 1) ]T exp 



n=l 



n 



sin nnr] 
sin 7r?7 



+exp - 




sin nr] J 



Let us note here the difference between the averaging procedure suggested above by 
relation (25) and the conventional scheme (see e.g. [[16]) for obtaining the characteristics 
of the electromagnetic radiation accompanying the motion of an ultra-relativistic particle 
in medium. The latter approach prescribes that, in the presence of multiple scattering, 
the radiative spectrum @ formed during the motion of a projectile in an external 
field along some particular trajectory must be averaged over all allowed trajectories. 
This operation is carried out with the help of the distribution function satisfying the 
corresponding kinetic equation [|l(|. Although the AIR radiation in medium is not an 
exclusion from this rule, in this case there is one peculiarity which justifies the use of 
the averaging procedure (25). As mentioned above, the motion of a positron channelled 
through a periodically bent crystal is characterized by two essentially different modes. 
The "fast" mode corresponds to the channeling oscillations inside the channel, the "slow" 
mode describes the motion of the particle along the channel centerline. The AIR is due 
to the "slow" mode and in the case a 3> d its characteristics are almost independent on 
the parameters of the fast channeling oscillations |2|, H. Hence, if one is interested in the 
spectral-angular distribution of the AIR only, then it is sufficient to totally disregard the 
channeling oscillations and to assume that the particle moves along the fixed trajectory 
which is defined by the profile of the bent channel centerline. The intensity of the AIR 
radiation is then defined by the number n of the A-periods which the particle passes 
before leaving the channel due to the multiple scattering. The value of n is subject to 
the probability distribution (|24"D which, in turn, is based on the exponential decay law 



Hence, to obtain the parameters of the AIR radiation in presence of the dechanneling 
it is sufficient to substitute the factor D N (rj) (see (|20|)) with its averaged value defined 
in (p7|). The sum on the r.h.s. of 
following expression for (Dn(t])): 

N\ fsmNnyy | 2iV 2 (1 + (2N d )~ L ) ^ _ ^ ^_AT 



(Dn(v)) 



exp 
N 



N,, 



is carried out straightforwardly resulting in the 

2iV 2 (1 + (2N,)- 1 ) 



smnr] 



1 + 4-/VJ sin 2 7irj 



(2f 



-exp — 



sin Nirrj 

NdJ 1 + 4A^| sin 2 7177 sin 7177 



2cos(A r + l)ixri + 



sin(A^ + 1)7777 
N d sin TiT] 



COS TIT] 
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Although more cumbersome, if compared with simple formula (|20|), the structure of 
r.h.s. of (28) can be readily analyzed. 

As well as the factor Dn{t]), the function (.0^(77)) has sharp maxima at the points 
rj — k — 1, 2, 3 . . . which correspond to the frequencies u>k defined as in ([H]). For given 
N the maximum value _Dat(/c), which is the same for all k, is equal to 

N 2 , for N/N d < 1 



N\ ( N 



(29) 

2N 2 , for N/Nd > 1 



and monotonously increases with the number of the undulator periods. 

The limit N/Nd <C 1 corresponds to a thin crystal, L <C L^. In this case the number 
of the channeling particles does not change noticeably on the scale L, so that the AIR 
intensity is proportional to N 2 as in an ideal undulator. For L ~ Ld the probability 
of the particle to channel through the whole crystal length decreases, and the value of 
{D^ik)) starts to deviate from the iV 2 -law. In the limit of very thick crystal, L 3> L^, 
{Dj^(k)) reaches its maximum value 2N% + Nd ~ 2Nj. 

The r.h.s. of (28) reduces to simple expressions in two limiting cases of thin and 
infinitely thick crystals: 



(D N (V)) 



'sin Ntttj^ 
sin7r?7 , 
2Nj 



for N < N d 

(30) 

for N > N d 



I 1 + ANj sin 2 -K7] ' 

Both limiting expressions have clear physical meaning. For iV -C Nd (L <ti Nd) the 
dechanneling is of no importance, and, therefore, the structure of the characteristic 
line of the AIR is that of an ideal undulator, eq. (|2Ti|). With iV (L) increasing the 
profile of the characteristic line changes and, finally, becomes of the Lorenz-type in the 
limit iV ^> Nd and for \r] — k\ 1. The qualitative explanation of this result can 
be given if one interprets the relation ( ^3|) in terms of quantum mechanics. Then, the 
r.h.s. is proportional to the squared modulus of the projectile wave-function which 
corresponds to the bound (channeled) state characterized by a complex energy. Indeed, 
the wave function of the channeled particle can be presented in the form (e.g. f24l ) 
^(i^t) = exp(— iet)ijj(r\\)if)(rj_, t), where ip(r\\) corresponds to the unbound longitudinal 
motion of the particle, and ^(r±,t) describes the transverse motion in the channel. 
Provided the multiple scattering is excluded the function ip(r±,t) = -^(rj_) exp(— ie^t) 
corresponds to a stationary bound state of the transverse motion (e± is the energy 
associated with the transverse degree of freedom). In this case, the normalization 
condition / dr_ L |'0(r, t ± )| 2 = 1 means, that the probability to find the particle inside 
the channel does not depend on time t, which, in turn, is related to the penetration 
distance through t = z/c. Random scattering by the electrons and nuclei can be 
incorporated in this picture by adding the imaginary term to the energy of the transverse 



14 



motion, e± — > e± + i(r/2)£, where T is the total width associated with the transitions 
to the unbound (dechanneled) continuum due to the multiple scattering. Hence, the 
normalization condition acquires the form J drj_|-0(r, t x ) | 2 = exp(— Ft) = exp(— Tz/c). 
Comparing this expression with (|23| ) one finds r = c/Ld oc 1/Nd which is exactly the 
quantity which defines the shape of the resonant line of (Dn(t))) in the case L 3> Ld 
(the second relation in (|30|)). 

The change in the shape of the characteristic line is illustrated in figure ^ where the 
function (D N {rf)) /Nj is plotted for several values of iV (as indicated) in the vicinity of 
the point 77 ~ k. For the sake of comparison we plotted also the dependence D^{rj)/N 2 
for N = Nd (the thin solid line) which characterises the line shape in an ideal undulator 
without the decrease in the volume density of the channeled particles. 



2.0 




N=N d /2 ~\ 


1.8 




N=N d j 

— - N=2N d j 


1.6 




N=3N d j 

N=20N d 


1.4 






1.2 


} 


\ ■: 
1 z 


1.0 


r / 


\ \ 

> ] 


0.8 


L , 

; / 


\ " 
\ : 


0.6 






0.4 






0.2 

A A 







-0.10 -0.05 0.00 0.05 0.10 



T]-k 

Figure 2. The dependences (Djv(r?)) /Nj (eq. (28)) for various N as indicated in 
vicinity of the point r\ = k, (k = 1, 2, . . .). The parameter Nd is set to 40. The thin 
solid line corresponds to D n (rf) / 'N 2 (see (^o|)) calculated for N = Nd- 

Hence, the expression (26) combined with (28) describes the spectral-angular of 
the total radiation in the range uo ~ uj}, uj c h formed during a passage of an ultra- 
relativistic positron along periodically bent crystallographic plane. All the complexity 
of the trajectory related to the channeling oscillations are incorporated in the undulator 
parameter p 2 (see (22)). The effect of the dechanneling is included in the factor (Djs[(rj)), 
eq. (28), which depends on two parameters, N and Ld. The number of the undulator 
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periods N can be varied by changing the length of the crystal L or/and by changing the 
period A of the shape function S(z) (see figure [I])). 

In the case of a straight channel the dechanneling length Ld depends on the type of a 
crystal and on the energy of projectile |15], |HJ . For a crystal bent with constant curvature 
1/R the quantity Ld becomes also dependent on R. The methods for estimating Ld for 



heavy projectiles channeling in crystals bent with 1/R = const are discussed in JT2 
To our knowledge, no detailed study has been carried out of the dependence of the 
dechanneling length on the parameters of the shape function S(z) in the case of a light 
projectile (a positron) channeling through periodically bent crystal. In the next two 
sections we present the approach and the results of numerical calculation of Ld in this 
case. 



4. Simulation of the dechanneling process in periodically bent crystals 

As it was demonstrated above the intensity of the radiation and the profile of the line 
of the AIR are defined, to a great extent, by the parameter Ld- 

This section is devoted to the description of the model which we used to calculate the 
dechanneling length Ld for an ultra-relativistic positron channeling through periodically 
bent crystal. Our approach is based on the simulation of the trajectories and the 
dechanneling process of an ultra-relativistic positron. This is done by solving the three- 
dimensional equations of motion which account for (i) the interplanar potential, (ii) the 
centrifugal potential due to the crystal bending, (iii) the radiative damping force, (iv) 
the stochastic force due to the random scattering of projectile by lattice electrons and 
nuclei. 

We analyze the dependence of L d on the energy of the projectile, the type of the 
crystal and crystallographic plane, the parameters of the shape function S(z). The 
latter include two factors: the parameter C and the ratio a/d (see eqs. (2-3)). 

The results of our calculations of the quantity L d , presented in section [|, are used 
further in section |6] to calculate the spectral-angular and the spectral distributions of 
the AIR. 



4-1. Equations of motion in periodically bent crystal with account for the radiation 
damping 

Below we outline the derivation of the equations of motion for an ultra-relativistic 
positron channeled in a periodically bent crystal. Contrary to the case of a heavy 
projectile, a light projectile (a positron, an electron) with e > 1600 m/Z looses its 
energy, when passing through matter, mainly due to the radiative losses (e.g. P3f ). 
In this inequality Z is the atomic number of the crystal atoms. Hence, the radiative 
losses exceed the losses due to the ionizing collisions starting with m 60 MeV in Si, 
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and ~ 11 MeV in W. Therefore, accurate treatment of the equations of motion of the 
channeled ultra-relativistic positron must account for the effect of the radiation damping 
IS , ^6], f|7| . In our scheme by introducing the classical radiative damping force . 

The equations of motion for a relativistic particle of mass m and charge e moving 
in an external static electric field E read as follows 



m^( 7 v) 

mc 2 7 



eE + f 

eE • v + f ■ v 



Here 



2e 3 
3mc 3 



7 (vV) E 



fvE) E - 



mc" 



7 



fvE) 



(31) 



(32) 



is the radiative damping force due to the presence of the electric field E. The 
latter is related to the interplanar potential U(r) through eE = —dU/dr, where 
r = y e y + z e z + x e x is the radius vector. 

Let the crystal centerline be bent as it is presented in figure |I[ We impose the 
following condition on the derivative of the periodic profile function S(z): 



S"{z) ~ (a/X) 2 = £ « 1 



(33) 



Then, the length L of the bent crystal, the interplanar spacing d, and the local curvature 
1/R(z) of the centerline satisfy the relations 

L = L (l + O(£ 2 )) « L , rf (l + O(£ 2 )) « do, JT 1 = |£"(2)|(l+0(£ 2 )) » |^(z)|(34) 

Here the subscript refers to the straight channel. When combined, the relations (H|) 
lead to the conclusion that the parameters of the bent channel, which are the channel 
width and the distance between the atoms of any crystallographic plane, are equal to 
the their values in the case of straight crystal provided one neglects the terms of the 
order £ 2 and higher. This, in turn, allows to formulate the following two conditions 
regarding to the interplanar potential in periodically bent crystal: 

1. within any single bent channel the potential U depends only on the variable 
p ~ y — S(z) which is the distance from the (x, y, z) point to the channel centerline 
(see figure ^) ; 

2. the dependence of U on p in the bent channel is identical to that in the linear one. 
Therefore, 

U(v) = U(p), p = y-S(z) (35) 
Eq. (|35D leads to the following definition of the interplanar static electric field E: 



eE = -U' 



S'e? 



(36) 



Using (36) in ( |3"2"D and in (31) and carrying out some algebra one ends up with the 
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Figure 3. The coordinate p in a bent channel: p is the distance from some point (y, z) 
in the channel to its centerline. To obtain p it is necessary to find the point (y , z Q ) on 
the centerline in which the tangent to the centerline is perpendicular to p. If omitting 
the terms S' 2 (z) and higher the p coordinate equals p — y — S(z). 



following system of coupled equations 

P= -(— + c*S» + ^pU») (37) 
yrrry m 3c j 

x = 77 (38) 

m7c z 

7= 2 \ U 39 

mc z y 3c m J 

Here r = e 2 /mc 2 = 2.818 x 10~ 13 cm is the classical electron radius. 

Since 1 — 7~ 2 = (x 2 + y 2 + z 2 ) / c 2 then dependence z(t) is uniquely defined provided 
the functions p(t), x(t) and j(t) are known. 

On the r.h.s. of (37-39) the terms of the order 7~ 2 are omitted. We also dropped 
out the terms of the order U/rwyc 2 . This is justified by the fact that for most of crystals 



the depth of the interplanar potential well Uq lies within 10. . .100 eV (see e.g. [16]). 
Therefore, for 1 ... 10 GeV positrons the parameter U^/m^c 2 ~ 10 -9 . . . 10~ 7 < 7 -2 . 
The profile of the periodically bent channel enters the equations (37-39) via the 
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function S"(z). Recalling (|34"D, the second term on the r.h.s. of (37) takes the form 
—c 2 /R(z) and, thus, represents the acceleration due to the centrifugal force. 

4-2. Equations of motion with account for the random collisions with target electrons 
and nuclei 

Random scattering of a projectile by target electrons and nuclei leads to deviation 
of the trajectory from that which is obtained by solving equations of motion (37-39). 
There are two different approaches aimed to account for these deviations. The first 
approach is based on the diffusion theory applied to describe the multiple scattering 



e.g. [|Tl|, [17|, H| |49fl . The second approach is in direct computer simulation of the 



scattering process |50], |12| . In the present paper we adopt the scheme similar to the one 



described by Biryukov [12 



Scattering from target electrons results in two changes in the projectile motion. The 
first one is the gradual decrease in the projectile energy due to the ionization losses. 
The second effect of the collisions is that they lead to a chaotic (random) change in the 
direction of the projectile motion. 

The scattering from target nuclei results in a chaotic (random) change in the 
direction of the projectile motion. 



4-2.1. Ionization losses for a channeling positron Although for an ultra-relativistic 
positron (7 ^> 1, v ~ c) passing through media the ionization losses are much smaller 
as compared with the radiative ones (see e.g. [[51], |52|) they are incorporated into the 
scheme which is described here. The mean energy loss per path ds ~ cdt due to the 
electronic scattering can be written as a function of the distance p from the midplane 

g, 0, [53 1 

, ~f\ / 2mc 2 T x 



r~ 1 — Kmc 2 {n e i) 

cdt . x 1 



ln- 



24 2 + ° {P) + M 



m- 



T 

J. rr 



(40) 



The notations used here £1X6 cLS follows: 

• K — 27irl k5x 1(T 25 cm 2 . 

• The quantity (n e i) is the mean concentration of electrons in the corresponding 
amorphous media. It is related to the mass density p m of the material, the 
atomic weight of the crystal A and the atomic number Z of its atoms through 
(n e i) = N \Z 'p m / 'A (with Na = 6.022 x 10 23 mol -1 being the Avogadro number). 

• n ei{p) is the local concentration of the electrons inside the channel. To calculate 
n e i(p) we use the Moliere approximation. 

• The quantity 5 is a so-called density effect correction |5_2], [S3 
relativistic positron one can use 5 ~ 21n7 



For an ultra- 
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The quantity T max is the maximum magnitude of the energy transfer from a 
projectile to a target electron. For an ultra- relativistic positron T max « 7771c 2 . 
The term C(p) is the correction due to the periodicity of the crystalline structure. 
The exact form of this correction was found in [ |55|j . 



Substituting e = •ymc 2 and T u 



7777c into 



we write 



-Kc (n a ) 



In 



mc 



I 24 (n e i) 



In 



7777c 



(41) 



To account for the influence of the ionization losses on the channeling motion one can 
add (|4lD to the r.h.s. of (39). Then, the resulting equation for 7 reads 



7 



mcr 



U'p + 



2r 7 2 V 




(42) 



4-2.2. Random change in the direction of motion due to ionizing collisions To account 
for the random change in the direction of motion due to single collisions with target 
electrons we followed the procedure proposed by Biryukov [|5(| for heavy projectiles 
but modifying his formalism for the case of a light projectile. For an ultra-relativistic 
positron travelling through a crystal the differential probability (per path ds = cdt) of 
the relative energy transfer p = (e — e')/e due to the the ionizing collisions with the 
quasi- free electrons is defined by the following expression |53] , 5B|: 

d/i ds 7 p l 

The relative energy transfer defines a round scattering angle 9, measured with 



respect to the instant velocity of the projectile, through (see e.g. f4"2 



cos 9 

P G [// 

7 



1 ~ /V Mm; 

1 -ap/p t 



mm; /^maxj > 
1 



i-i 

7 



a 



7 + 1 



i-i + o( 7 -) 



(44 a) 

(446) 
(44 c) 



Here J ~ 16 Z 0,9 eV characterises the average ionization potential of the crystal 
atoms. The maximum relative energy transfer p = /i max results in 9 = n/2, whereas 
P = Pram *C 1 leads to nearly forward scattering. 

Following Biryukov |5(J we used the distribution (f£|) to generate (randomly) the 
magnitude of 9 when integrating the system of equations (37-38) and (42). The 
magnitude of the second (the azimuthal) scattering angle is not restricted by any 
kinematic relations, and is obtained by random shooting (with a uniform distribution) 
into the interval [0, 2ir]. 
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The algorithm of random generation of the 0-values is as follows. The two-fold 
probability dP/ d/z ds satisfies the normalization condition 

/ ^ion 

/ / T^r d ' uds = 1 ( 45 ) 

JO Jfimin dfl dS 

where L; on is the interval inside which the probability of the projectile to undergo the 
ionizing collision accompanied by arbitrary energy transfer is equal to 1. The parameter 
L ion (defined from (|5|)) is expressed through /i m i n , /i ma x and the local electron density 
n e i(p) as follows 

hon = > (46) 

/ A^min 

Using these notations eq. (f|^) can be written as 

d p = p- W{ii) dp W{ii) = ^ (47) 

Aon A* 

Here the factor As/L; on defines the probability of the collision (any) to happen on 
the scale As, whereas the factor Win) d/i represents the normalized probability of the 
energy transfer between \i and \i + d/i. To generate random deviate with the probability 



distribution Win) we used the algorithm described in []57| . 

The scheme outlined above implies that the ionizing collisions are treated as events. 
Hence, for each step As = cAt of integration of the the system (37-38, 42) we first 
simulate the probability of the event to happen by generating a uniform random deviate 
x G [0, 1] and comparing it with As/L ion . If x < As/L ion then the scattering angles 9 
and <p are calculated and used to modify the direction of motion of the projectile but 
leaving the magnitude of the projectile velocity unchanged. 

4-3. Random change in the direction of motion due to scattering from nuclei 

The change in the projectile direction of motion due to the collisions with crystal nuclei 
was taken into account at every step when integrating the system (37-38), (42). The 
square of small deflection angle A8 2 per path As, was computed from a Gaussian 
distribution 

— = ; exp = , k — — — (48) 

Here k 2 = 9 2 / As is the mean-square deflection angle due to the scattering from nuclei 



in crystals. In the approximation suggested by Kitagawa and Ohtsuki [58 



where (n n ) is the mean concentration of nuclei in an amorphous media, n n (p) is the 
local concentration of the nuclei inside the channel. With the thermal vibrations taken 
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into account the quantity n n (p), corresponding to the distribution of the nuclei of two 
neighbouring planes versus the distance from the midplane, is given by 



n n (p) d 



(d/2 + p) 2 \ , / (d/2-pf 



6XP l j+^V 2u% 



(50) 



M ^2nu 2 T 

with T standing for the crystal temperature and ut denoting the thermal vibration 
root-mean-square amplitude. 

The quantity (k 2 ) on the r.h.s. of (|49|) is the mean-square of A8 2 /As in an 
amorphous media. For an ultra- relativistic projectile it is given by [|17| 

<* 2 > = ~§ ~T~ (51) 



e 2 L 



rail 



where e s = 21 MeV, and L ra( i is the radiation length |5I 



« 4« r 2 Z 2 (n n ) In (l83 Z^ 3 ) (52) 
5. Numerical results of the dechanneling process simulations 



We used the scheme outlined in sections fO] and JO] for computer modelling of 
the dechanneling process of ultra-relativistic positrons in periodically bent crystals. 
The goal of these calculations was to obtain realistic dependences of the number of 
channeled particles n(z) versus the penetration distance z. The interplanar potential 
was considered within the the Moliere approximation at the temperature T = 150 K. 

The calculations were performed for 5 GeV positrons channelling along the (110) 
crystallographic planes in Si, Ge, and W crystals. The shape function S(z) was chosen 
in the form (|TED. In this case the minimum curvature radius of the centerline bending 
equals -R m i n = l/(ak 2 ) = \ 2 /47i 2 a, and, consequently, the parameter C defined in (2) 
reads as C = eak 2 /U^ ax . Thus, for given crystal and crystallographic plane, and for 
given e the parameter C and the ratio a/d uniquely define the period A. 

In our calculations C was varied from to 0.5, and the considered ratios a/d are: 
(the case of a straight channel), 5,10, 15 and 20. 

For each pair of the parameters C and a/d we simulated 2000 trajectories by solving 
the coupled system (37-38) and (42), which was integrated with the account for the 



random scattering from target electrons and nuclei as explained in section |4.2| . The 
initial conditions were chosen as follows. The coordinates x and z were set to zero 
at the entrance, the transverse coordinate p was obtained by random shooting (with 
a uniform distribution) into the interval [— d/2, d/2]. The vector of initial velocity vo 
was aligned with the tangent to the centerline at the entrance but with allowance for 
the spread in the incident angles 9 G [— 6 max , 6 max ]. Hence, the initial values of the 
velocities x/c and p/c were obtained by random shooting (with a uniform distribution) 
into the interval [— O max , ©max]- For each crystal the parameter O max was chosen as 
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10 4 0l- The values for the Lindhard angle 0l for 5 GeV positron channeling along 
(110) planes in Si, Ge and W are equal to 9.57 x 10~ 5 ,1.28 x 10~ 4 ,2.35 x 10" 4 , rad, 
respectively. Thus, the considered spread in the incident angles is about 10~ 8 rad which 
is achievable in modern high-energy e + e~ colliders (see e.g. |5"6fl). 

The results of calculations are presented in figures ||-|| [5] and in table IB. 

The first three figures represent the dependences of n(z) (normalized to the entrance 
value n(0) = n ) versus the penetration distance z for Si, Ge and W crystals. 
The ratio a/d equals 10. The curves in the figures refer to different values of C 
as indicated. The value C = stands for the case of a straight channel. The 
corresponding values of the spatial period A of the shape function S(z) can be calculated 
as A = 2tt\J (ed/U{ QSLX )(a/d)C~ 1 and are presented in the table. 
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Figure 4. The calculated dependences n(z)/n(0) versus penetration distance z for 
5 GeV positrons channeling along the (110) in Si crystal for various values of the 
parameter C (see (2)) as indicated. The data correspond to the shape function 
S(z) = asin(27rz/A). The a/d ratio equals 10. For each indicated C the corresponding 
values of A, and the calculated magnitudes of the dechanneling lengths L c d and the 
number of undulator periods N% = L c d /\ are presented in Table |l|. See also the 
commentaries in the text. 



Let us discuss these results. First to be mentioned is that in all cases the dependences 
n(z)/n(0) are monotonously decreasing functions which are smooth for all z except for 



23 



1.0 f 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
penetration distance z, cm 

Figure 5. Same as in fig. ||| but for Ge crystal. 



the steep change in the magnitude in the region close to z — 0. This change, An , 
is the smallest for the straight channels and its magnitude increases with C. Such a 
behaviour is explained as follows. The probability for the projectile to undergo the large- 
angle scattering from the crystal electrons or/and nuclei is proportional to the volume 
densities n e i(p) and n n (p) (see ( ^6H50| )). The large-angle scattering, i.e. when 6 > 0l, 
results in the increase of the kinetic transverse energy, which becomes high enough for 
the projectile positron to leave the potential well. Hence, all those particles which at the 
entrance move in the region of high nuclear and electron densities leave the channeling 
mode almost immediately. The volume density of crystal electrons is comparatively 
small in the inner region of the channel and sharply increases as \p\ > d/2 — a^p (here 
otf = 0.8853Z _1 / 3 ao is the Thomas-Fermi radius of the crystal atoms, a is the Bohr 
radius), the density of target nuclei becomes noticeable for \p\ > d/2 — ut (mt standing 
for the thermal vibration root-mean-square amplitude). This is illustrated by figure 
[7|. Since in our simulations we assumed that the particles are uniformly distributed in 
the coordinate p at the entrance point, then the fraction ~ n(0) (d — 2a^Y)/d leaves 
the channel close to z = 0. This the mechanism of dechanneling in the vicinity of the 
entrance point [|12[] is intrinsic for both straight (C = 0) and bent (C > 0) channels. 
In addition to it, in the case of bent channel there is a process which is called bending 
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Figure 6. Same as in fig. ||| but for W crystal. 



dechanneling ||12f , and which is due to the change in the depth of the potential well. 
Indeed, inside the bent channel the transverse motion of the particle is subject to the 
action of the effective potential (see (37)) 

U eff (p) = U(p) - eS"{z) p = U(p) - p (53) 

where and R(z) is the (local) curvature radius of the channel. For a channel bent as 
described by ( |16D R(z) = [(\/2n) 2 / a] sin(27rz/A). The particle could be trapped into 
the channelling mode provided its total energy associated with the transverse motion is 
less the minimal value, t^ e //(Pmin), of the two maxima points of the asymmetric potential 
well described by (^3]) |7], 0. The potential U e jf(p) reaches the magnitude of U e ff(p min ) 
at some point p m i n which satisfies the condition |p m i n | < d/2 and the absolute value of 
Pmm decreases with the growth of a. Hence, the more the channel is bent, the lower 
the allowed values of the channelling oscillations amplitude are, and, consequently, the 
narrower is the interval of the initial po values for which the trajectory from the very 
beginning will correspond to the channeling one. 

Another feature to be mentioned in connection with the curves in figures |]-[)] is 
the noticeable decrease, with the growth of C, of the penetration distance which the 
projectile can achieve. If one assumes that the dependence n(z) satisfies the exponential- 




Figure 7. Volume densities of electrons and nuclei versus the distance p from the 
midplane calculated within the Moliere approximation (at T = 150 K) for (110) 
channels in Si, Ge and W crystals. The dashed lines mark the ratios (d/2—a,TF)/d. The 
values of d and a TF for Si, Ge, W are equal 1.92, 2.00, 2.45 A and 0.194, 0.148, 0.112 
A, respectively. 



decay law (see (|23|) ), then the curves demonstrate that the dechanneling length is 
the decreasing function of C. This effect has been discuss in literature ( Hl2|, |i"3fl ) 
in connection with the heavy particles dechanneling in crystals bent with constant 
curvature. Physically clear quantitative arguments based on the relationship 

between the diffusion coefficient and the mean-square angle of multiple scattering lead 
to the following estimate for the dechanneling length [12, [| in the channel bent with 
the mean curvature 1/R 



U&R) = (l - RJR 



-\2 



L e d (oc) 



L e d (oo) 



256 otf d 



9tt 2 r 



E 7 



(54) 



here the quantity R c = s/U^ x is the critical (minimal) radius consistent with the 
channelling condition in a bent crystal (2). The quantity L c = In ^y^ymc 2 //) —23/24, 
with / = 16Z 0,9 eV standing for the Thomas-Fermi ionization potential of the crystal 
atoms, is the Coulomb logarithm characterizing the ionization losses of an ultra- 
relativistic positron in amorphous media with account for the density effect (see e.g. 

00). 
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The quantity L e d (co) is the estimated value of the dechanneling length for a straight 
channel. The factor (l — R c /Rj is the correction due to the decrease of the potential 
well of U e ff(p) (see (|53|)) in bent crystal. 

In periodically bent crystal the mean curvature radius is proportional to -R m i n (in 
particular, for the shape (|16D, R = 7ri? min /2), hence the ratio R c /R oc C. This estimate 
was used in [^|, [7], |£| to calculate the characteristics of the spontaneous and stimulated 
emissions based on the AIR effect. 



Basing on the estimate (|54[) one deduces that the dechanneling length is mainly 
determined by the value of the parameter C and is not influenced by the a/d ratio. This 
conclusion is supported by the results of our calculations. In figure (§) the dependences 
n(z)/n(0), obtained for different C and a/d values, are compared. It is seen that the 
variations between the curves with different a/d but the same C is much smaller than 
for the curves with the same a/d ratio but different C values. 
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Figure 8. The calculated dependences n(z)/n(0) versus penetration distance z for 
5 GeV positrons channeling along the (110) in Si crystal for various values of the 
parameter C and of the ratio a/d as indicated. The data correspond to the shape 
function S(z) = asin(27rz/A). 



The results of our calculations of the dechanneling lengths for 5 GeV positrons in 
Si, Ge and W crystals obtained at different values of the parameter C but for fixed 
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ratio a/d = 10 are presented in table [I]. They are notated through L c d (the quantities 
N d = L c d /\ are the number of the undulator periods on the scale of L c d ). The presented 
values of L c d were obtained by fitting the smooth part of the curves n(z)/n(0) (i.e. for 
the distances z beyond the region of the steep decrease) with the exponent exp (—z/L d ). 

For each C we included also the corresponding values of the undulator period A, the 
minimum curvature radius -R m i n = 1/ (ak 2 ) = A 2 /47r 2 a, the energy of the first harmonics 
of AIR, huii, obtained from (^TJ) for the forward emission ($ = 0), and the values of the 
undulator parameter p. 

The quantities L e d represent the the dechanneling length estimated from (|54]) , 
N d = L e d / A is the corresponding number of the undulator periods on the scale of L e d . 

Comparing the calculated values L c d with those estimated from ( |54l) one notices that 
the latter underestimates the dechanneling length in the case of small values of the 
parameter C: C < 0.25 for Si, < 0.20 for Ge, and < 0.15 for W. It happens mainly 
due to the fact that the quantity L d (oo), which enters (|54D, was obtained by using the 
Lindhard approximation for the interplanar potential U(p). The latter overestimates 
the electron density inside the channel as compared with the Moliere approximation 
used in the present calculations. Thus, the average multiple scattering angle is higher 
if calculated within the Lindhard approximation, and, consequently, the dechanneling 
lengths are lower. With C increasing the ratio L c d / L e d becomes less than one. This is 
explained as follows. When the parameter C, and, consequently, the curvature of the 
channel bending, increases, then the minimum of the effective potential ([53]) is shifted 
from the midplane towards the atomic plane. Hence, the channeled particle moves in 
the region of higher electron density and, therefore, the probability of the large-angle 
scattering increases ]T2| , |T3] ] . This is mechanism of the dechanneling, which is additional 
to the bending dechanneling discussed above, is not accounted for by the estimate fl54]) . 
Comparing the L c d values for various C it is seen that for large C the dechanneling length 
rapidly goes to zero, and the main reason for this behaviour is the increase of the role 
of the large-angle scattering. 

Another important feature to be mentioned is that the estimate (|54j) predicts the 
decrease in the dechanneling length for heavier crystals, whereas the more accurate 
calculations within the framework of the Moliere approximation demonstrate that for 
C G [0 . . .0.1] the largest dechanneling lengths are in W, followed by the germanium 
crystal, and finally, the lowest values of L c d are in the case of channeling in Si. With 
C increasing the situations is reversed. To explain this behaviour let us compare the 
cases of Si and W crystals. The maximum values of the interplanar potential U(p) at 
temperature T = 150 K are equal to 22.9 eV for Si, and to 138.6 eV for W. Thus, the 
critical transverse energy of the channeled positrons is approximately 6 times higher in 
W. Hence, were the mean electron densities in these crystal equal then the dechanneling 
length in W would be 6 times higher as compared to the silicon case. The ratio of 
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Table 1. Dcchanncling lengths for 5 GeV positron channeling along the (110) planes 
for various crystals and for various values of the parameter C (eq. (2)). The data 
correspond to the shape function S(z) = a sin(27rz/A). The a/d ratio equals 10 except 
for the case C = (the straight channel) . The quantity L c d presents the results of our 
calculations, N d = L c d /\ is the corresponding number of the undulator periods, L e d is 
the dechanneling length estimated according to (|54|), N d = L d /X. Other parameters 
are: d is the interplanar spacing, R c = e/C^ ax is the critical (minimal) radius consistent 
with the condition (2), Tiu>i is the energy of the first harmonic of the AIR for the forward 
emission (see (21)), p is the undulator parameter,. 
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the mean electron densities for the inner regions of the (110) channels in W and Si 
is approximately 3.8. Therefore, the realistic estimate for the ratio (L d ) w / (L d ) Si is 
6/3.8 ~ 1.6 which is very close to that which can be obtained from table [I] for low 
C-values. Our results for the relative magnitudes of L d for a positron channeling in 
Si and W coincides with the conclusion of made in [^] where the calculations of the 
dechanneling lengths were carried out for heavy particles in the case of moderate bending 
of the crystals. 

As C increases, and the trajectory of the particle shifts towards the atomic plane, 
the electron density in the region of the projectile motion increases as well. This increase 
is more pronounced for W (see figure (|7|)) leading to the relation (L d ) w / (L d ) Si < 1 in 
the range C > 0.15. 

The important conclusion which can be made on the basis of the results presented in 
this section is the following. To achieve higher yield of the monochromatic radiation by 
means of the crystalline undulator it is necessary to use heavier crystals and to restrict 
the parameter C by the range 0.05 . . . 0.15. Indeed, for these C values the number of 
the undulator periods iV « N d ^> 1 (this is valid for all crystals presented in the table), 
but the intensity of the AIR, which is mainly defined by the factor N 2 Aj (see (26) 
and (|30|) ), will be almost by the order of magnitude lager in the case of W than for Si 
and/or Ge crystals. 

6. Results of calculations of the spectral-angular and spectral 
distributions of the AIR 

The simulation of the dechanneling process, described in the preceding section, allows 
to establish the ranges of the parameter N d which to a great extent defines the 
characteristics of the AIR radiation in the presence of the dechanneling. In this section 
we present the results of calculations of the spectral-angular (see eq. (26)) and the 
spectral distributions of the AIR. 

The general features of the angular distribution of the radiation formed in an 
acoustically based undulator were discussed in where the effect of the dechanneling 
was taken into account by restricting the length of the crystal (and, hence, the number 
of the undulator periods) by the value L d , which, as it was demonstrated above, in 
many cases underestimates the dechanneling length, and, consequently, leads to the 
underestimation of the photon yield. 

Figures |](a)-|l(d) present the spectral distribution (26) of the radiation emitted 
along the undulator axis, h~ 1 (dE N /duj dr2 n )^ =0 °, for 5 GeV positron channeling along 
(110) planes in Si and W crystals. The spectra were calculated by using the parameters 
A and p as indicated in table [I] in the lines C = 0.15 for Si, and C = 0.05 for W. The 
number of the undulator periods varied from N = 4A^, as in figures |5|(a,c) and the 
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solid curves in figures |9](b,d), to N = N%/2 (the long-dashed lines in figures [J](b,d)). 
For fixed A different values of iV can be achieved by changing the length of the crystal. 
The sinusoidal profile of the channel bending can be created by applying the transverse 
acoustic wave. The velocities of the transverse supersonic waves propagating along (110) 
directions are equal 5.84 cm/s for Si [50| and 5.17 cm/s for W crystal pl |. Thus, the 
corresponding frequencies of the acoustic wave needed to to obtain A = 63.0 /xm in Si 
and A = 42.2 /iin in W (see the table) are: u Si = 93 MHz, z/ w = 123 MHz. 
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Figure 9. Spectral distribution (26) (in 10 9 sr _1 ) at i9 = for 5 GeV positron 
channeling along periodically bent (110) planes in Si (figures (a) and (b)) and W 
(figures (c) and (d)) crystals. The parameters of the crystals bending are as presented 
in the table for C = 0.15 (for Si) and C = 0.05 (for W). The upper figures (a) and (c) 
reproduce (di?/v/?id(jj d£ln) in the wide ranges of Li and correspond to N — 4iV^. The 
numbers enumerate the harmonics (in the case of the forward emission the radiation 
occurs only in odd harmonics). The profiles of the first harmonic peak (figures (b) and 
(d)) are plotted for N = 4N% (solid lines), N = 2N% (dotted lines), N = N% (dashed 
lines), N = (long-dashed lines). In both figures the thin solid line corresponds to 

the non-averaged spectrum (jT^) calculated with the number undulator periods N = iVJ 
which follows from (p4). 



The upper figures illustrate the structure of the spectral distribution in wide ranges 
of the emitted photon energy. Each peak corresponds to the emission into the odd 
harmonics [[24], |TB|, the energies of which follow from ([H]) if putting $ = 0°. It is seen 
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that all harmonics are well separated: the distance 2hu>i between two neighbouring peaks 
is 2.74 MeV for Si and 1.78 MeV in the case of W, whilst the width of each peak HAu, 
estimated from (0) (the case N > iVJ), is « 8.7 keV for Si and « 2.5 keV for W. The 
intensity of the first-harmonic peak in W is approximately 4.5 times larger than in the 
case of silicon crystal. This ratio is in agreement with the expression (26). Indeed, the 
intensity of the first-harmonic peak is proportional to the factor p 2 (ui/uo) 2 (Dn(1)), 
which, in turn, is proportional to p 2 /(2 + p 2 ) 2 (iVj) 2 (see eqs. fl2"T| ) and (PD|)). Using the 
data presented in the table one obtains that the ratio of these factors calculated for W 
and Si is ~ 5. 

The difference in the magnitudes of the undulator parameters, p = 1.87 for Si and 
p = 3.26 for W, explains number of the harmonics visible in the spectra [p4| , £L£|. 

Figures ||(b,d) exhibit, in more detail, the structure of the first-harmonic peaks. For 
the sake of comparison we plotted the curves corresponding to different values of the 
undulator periods. It is seen that for N > N% the intensity of the peaks is no longer 
proportional to N 2 , as it is in the case of the ideal undulator without the dechanneling 
of the particles (Section ||). For both Si and W crystals, the intensities of the radiation 
calculated at iV — > oo exceeds those at N = 4A^ (the thick solid curves in the 
figures) only by several per cent. Thus, the solid lines correspond to almost saturated 
intensities which are the maximal ones for the used crystals, projectile energies and 
the parameters of the crystalline undulator. It is worth noting that these intensities, 
obtained with realistic values of the dechanneling lengths L% are noticeably higher than 
the ones corresponding to iV > iVJ (the thin solid lines in the figures). The latter were 
used in for the estimations of the yield of the AIR, both spontaneous and stimulated. 
It it seen from the figure that our earlier estimations were adequate for light crystals 
(the maximum of the thick solid curve exceeds that of the thin one by a factor of 3.5) 
but they essentially underestimated the yield of the AIR in heavy crystals, where the 
ratio of the intensities is almost 10. 

In figures |iy and we present the results of the spectral distribution of the radiation 
h^ 1 (dE^/duj) obtained by integration of the r.h.s. of (26) over the emission angle dfi n . 
The calculations were performed for 5 GeV positrons channeled in Si and W with the 
bending parameters as those mentioned above. 

It is known from general theory of the undulator radiation by an ultra-relativistic 
projectile |24], [16] that it is emitted into the narrow cone along the undulator axis. 
The opening angle of this cone, i9 c depends on the two parameters, namely, on the 
relativistic factor 7 and on the undulator parameter p. For p 2 <C 1 (the so-called dipole 
case) the radiation is concentrated within the cone with the opening angle ~ I/7. In the 
opposite limit, p 2 ^> 1 (the non-dipole case) the opening angle becomes larger, ~ v/l- 
Uniting both limiting cases one can state that for given 7 and p the maximum angle of 
the emission can be estimated as d c ~ max {7~ 1 ,P7 -1 } <C 1- In the case of a planar 



32 



0.8 

0.6 

0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 



- 1 1 1 


i i i i i i - 

i 

1 3 

1 , ', , I , ; , l 1 




(b) : 

f\ /I [\ rA } 


1 I 1 ' 

: 


(c) J 

f\ rA rM i — 


1 1 

/ 


(d) : 


i i i i i i i i i i f, 

; / 


(e) : 



0.0 1.0 2.0 3.0 4.0 5.0 

photon energy, MeV 



6.0 



7.0 



8.0 



Figure 10. Spectral distribution (55) for 5 GeV positron channeling along periodically 
bent (110) planes in Si. The five graphs correspond to different values of t9 max : 
(a) tf max = 0.10c, (b) tf max = 0.2 C , (c) z9 max = 0.3 C , (d) tf max = 0.4 C , (e) 
= 0.5 $ c , with $ c = p/"f — 0.187 mrad. The numbers in figure (a) enumerate the 



harmonics. 



undulator, which is considered here, and for p > 1 the cone is asymmetric with respect 
to the azimuthal (<p) angle of the photon emission 0. The opening angle is the largest 
for the emission within the undulator plane (<p = 0°) and reaches its minimum value for 
the emission into the perpendicular plane (<p = 90°). 

The five graphs in each of the figures reproduce the dependence of the quantity 

on io for different values of the parameter i? max as it is indicated in the figures. 

The parameter i? c , defined above, equals to 0.187 mrad for Si and 0.326 mrad for 
W crystals. 

The graphs (a)-(e) in each figure illustrate how the pattern of the spectral 
distribution changes with enlarging the opening cone tf max . For small "# ma x, when only 
the radiation emitted in the nearly forward direction is taken into account (the graphs 
(a) in the figures), the shape of the distribution (55) is close to that exhibited in figures 
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Figure 11. Same as in figure ( |ic| ) but for tungsten crystal. The parameter i3 c = vll 
equals to 0.187 mrad. 



|9] (a) and (c). The radiation into the odd harmonics dominates over the even-harmonics 
peaks which, nevertheless, become visible in contrast with the case of emission at — 0. 

As the cone of emission becomes larger (the graphs (c)-(e) in the figures) the width 
of the peaks grows and their shape becomes asymmetric. The enhancement of the 
width follows from eq. (plf) which connects the harmonics frequencies with the emission 
angle, and from figure ^|. For a given number k of the harmonic its frequency 
is the decreasing function of the emission angle, so that the center of the peak shifts 
towards lower values of the frequency as $ grows. Hence, the integration of (26) over 
the interval 1?, 1? + A$ leads to the appearing of the emission within the frequency range 
~ [a;fc(i?) — Au, cufc(^)] where 

4 7 2 $ a?? 



2 + p 2 + 2 7 2 ^' 



We illustrate this estimate by figure where the graphs (a)-(e) correspond to the 
contributions of different intervals of the emission angle to the spectral distribution 
of the AIR for 5 GeV positron in tungsten. It is clearly seen that each interval of 9 
corresponds to the emission radiated within particular ranges of u>. 
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Figure 12. Spectral distribution of the AIR for different intervals of the emission 
angle ■&: (a) < d < 0.1 C , (b) 0.1 i? c < d < 0.2 i? c , (c) 0.2 i? c < i? < 0.3 i9 c , (d) 
0.3 i? c < 1? < 0.4 (e) 0.4 < i9 < 0.5 t9 c , with i? c = 0.187 mrad. The graphs 
correspond to 5 GeV positron channeling along periodically bent (110) planes in W. 
The numbers in figure (a) enumerate the harmonics. 



In figures [10| and [11] the graphs (e) correspond to the value $ max = 0.5 ?9 C for the 
upper limit of integration over the emission angles. It is seen that initially well-separated 
narrow peaks (the graphs (a)) have merged, except for the first three harmonics. The 
further increase of i? miffi up to ^ c produces the continuous spectrum of radiation although 
with the irregularities in the vicinities of Uk(0). 



7. Conclusions 



In this work we have described the general formalism and the algorithm for effective 
calculations of the characteristics of the AIR as well as of the dechanneling length of 
positrons in periodically bent crystals. 

The calculations performed in the present paper confirm the main result of our 
previous considerations concerning the possibility to create a powerful and easily tunable 
source of monochromatic radiation in the X or/and 7 range by means of crystalline 
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undulator. They also demonstrate that the knowledge of the accurate values of the 
dechanneling lengths for various projectile energies, crystals and crystallographic palnes, 
and various parameters of the profile function S(z) is essential for obtaining reliable data 
on the intensity of the AIR. 

The numeric results presented above refer to a particular shape of the bent channel 
which can be achieved by applying monochromatic transverse supersonic wave. However, 
both the formalism and the computer code allows to investigate, by comparatively simple 
means, the radiative spectra formed during the channeling of ultra-relativistic positrons 
along the crystallographic planes whose shape is described by arbitrary periodic shape 
function. 

In this connection we want to mention the problem concerning the possibility to 
generate the AIR photons not in the range of up to several MeV but much higher. It can 
be achieved by using positron beams of higher energies, e — 10 . . . 100 GeV. The increase 
of e results in the increase of both the emitted photon energies and the dechanneling 
length. The latter leads to the increase in the intensity of the AIR. The obvious difficulty 
of the stable work of the crystalline undulator in this range of e is in sharp increase of 
the radiative energy losses which are mainly due the ordinary channeling radiation. The 
radiative losses lead to a noticeable decrease of the projectile energy during the passage 
through the undulator. Hence, the use of the shape function with constant parameters 
a and A does not result in the intensive monochromatic undulator radiation. However, 
we believe, that this problem can be overcome by considering S(z) with the varied 
parameters, a(z) and X(z). Then, provided the dependence e(z) (z = ct) is known it 
is possible to adjust the parameters a(z), X(z) so that the characteristic frequencies of 
the AIR will not change as the projectile penetrates through the crystal, and, therefore, 
the monochromaticity will be restored. This work is in progress at the moment and will 
become the subject of another publication in the near future. 
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